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Abstract 

With the diquark structure of A c , we investigate the branching 
ratios of A c — > nii + and pir . The results show that without 
considering the final state interaction (FSI), the branching 
ration of A c — > pn + is only of order 1CU 6 whereas taking into 
account the FSI effects, this ratio could reach 10~ 4 and is 
at the same order as A c — > nn . Concrete values depend 
on phenomenological parameters adopted in the calculations. 
These branching ratios can be measured in the experiments 
to come. 



PACS numbers: 14.20.Lq, 13.30.-a, 12.39.Hg, 12.39.-x 



I. Introduction 



In the recent years, there are remarkable progresses in heavy flavor physics for both 
theoretical framework and applications. As long as heavy flavors are involved in the 
concerned physical processes, due to an extra symmetry SU(2)f <g> SU(2) S , calculations 
become simpler and more reliable, and the physics picture is clearer as well ffl. Based 
on the understanding, the heavy quark effective theory (HQET) has been established 
and applied to investigate the physical processes. In the HQET, at the leading order of 
1/rriQ expansion, the non-perturbative QCD effects are attributed into one form factor, 
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i.e. the Isgur-Wise function. With it, one can make predictions and extract significant 
information from data, such as the CKM entries. 

Heavy-meson related processes have attracted intensive interests because a great amount 
of data about B and D have been accumulated so far. By contrary, the processes of heavy 
baryons have not received as much attention as for mesons. The reason is two- fold. First, 
there are not sufficient amount of data on heavy baryons available yet. Secondly, one has 
to deal with a three-body problem theoretically which is much more complicated than 
the meson case. On the other side, there are many unsolved problems for baryons, for 
example, the lifetime of A& is 27% shorter than that of B-meson, and it indicates the 
shortcoming of the spectator scenario ||. 

The diquark picture was proposed a long time ago and it has been applied to study 
many processes involving baryons [|J. Obviously, for a baryon containing two heavy and 
one light quarks, the two heavy quarks constitute a diquark which serves as a color-field 
source for the light flavor ||. Instead, for the one- heavy and two-light quark structure, it 
is suggested that the two light quarks can also form a light diquark. 

In this work we are going to employ the diqaurk picture to study the Cabibbo- 
suppressed A c decays. The reason for doing such calculations is that in the near future 
there will be a great amount of data on A c decays available and probably, higher statistics 
for the rare decay modes will also be achieved, so we can use the data to testify our theory. 
Moreover, in these decays the final state interaction (FSI), which is also a very important 
subject, is more important than that in the case of In this work, we will investigate 
how much the FSI effects could change the branching ratios. 

We will mainly use the factorization ansatz which is not accurate for meson decays, 
as many authors pointed out || |^, |8|]. Even though the non-factorizable effects which 
are proportional to the color-octet currents may bring up some errors for the theoretical 
predictions, later we will show that they do not affect our qualitative conclusion. Besides, 
there are some uncertainties coming from other approximations and all of them are due 
to lack of sufficient knowledge on the non-perturbative QCD. 

All the weak effective vertices employed in this work are standard |§. Then we can 
calculate the weak transition in a straightforward way up to the hadronic matrix elements. 
To evaluate the hadronic matrix elements, one needs to invoke models, because they are 
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fully governed by the non-perturbative QCD. Since we are dealing with transitions from 
A c to lighter baryons A,p,n, etc., we cannot simply use the Isgur-Wise function which 
applies only to the transitions between heavy flavors. In terms of the diqaurk picture, 
we employ the wave-functions of baryons which are presented in literature to carry out 
calculations of the corresponding hadronic matrix elements. They finally are attributed 
into two independent form factors. 

Another important issue is the FSI effects which may bring up dramatic changes to 
the predicted branching ratios. As well known in the meson case, FSI plays a remarkable 



role and sometimes is almost crucial fL0| . Thus one can believe that in the baryon case, 
FSI may also be significant. In fact, our numerical results show that without considering 
FSI, the branching ratio of A c — > prr is only of order 10~ 6 . However, A c — > pir can be 
realized via A c — > AK(K*) — > pir or A c — > nir + (p + ) — > pre inelastic final scattering, and 
the resultant branching ratio could reach order of 10~ 4 . The results can be tested in the 
future experiments. 

Obviously our calculations depend on some phenomenological ansatz and certain pa- 
rameters, but all what we adopt in this work have been widely used in literature and most 
of them have been tested by experiments to some extent, hence they are somewhat trust- 
worthy or at least reasonable. Definitely, our results depend on the chosen parameters, 
but not very sensitively. In other words we do not fine-tune the numbers, therefore one 
can trust the order of magnitude of the results. 

The work is organized as follows. After the introduction, we present our formulation 
and meanwhile very briefly review some necessary knowledge. In Sec. Ill, we give our 
numerical results as well as the adopted parameters and the last section is devoted to our 
conclusion and discussion. 

II. Formulation 

The weak effective hamiltonian is written as (for c — ► q transition) || O] 

^// = %E^MC i (/i)o J , (1) 

V 2 i=i 

where q (i = 1, 6) are the QCD modified Wilson coefficients and the operators 0« have 
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the following expressions: 

Ox = Ua-f^l - 75)g/3g/37 M (l - 75)C Q , 2 = UJ^l ~ 7fe)?ffT M (l ~ ls)c, 

3 = u 7m (1 - 75 )c^gV(l -7s)g', 4 = u a j„(l - j B )c p £ - 75 )g^ 

<?' <?' 

5 = «7 M (1 -75)c^g'7 M (l +75)9', C 6 = m q 7 m (1 -75)^^ g^7 M (l + 7 5 )g^, 

(2) 

where a and /? are color indices, and q' = u, d, s. In the Hamiltonian we have omitted 
the operators associated with electroweak penguin diagrams. 
The corresponding transition matrix elements are 

< f\H eff \i >= %Y.Vh KM c^) <m\i>. (3) 



In the limit m c — ► 00, the baryonic transition matrix element for A c — ► / can be formed 

as 

< /(p/)|?7m(1 - l 5 )c\A c (v) >= U f {p f )[F x {v ■ p f ) + fF 2 {v ■ P/ )] 7/1 (l - 75)^00, (4) 

where t> is the four-velocity of A c , q = u,d, s, and Uf and U\ c are Dirac spinors for / and 
A c , respectively. The two form factors F\ and F 2 must be evaluated in concrete models. 

To calculate the transition amplitude, it is usually assumed that the process only oc- 
curs at the heavy quark arm while the light flavor quark (diqaurk) stands as a spectator. 
Recent research points out that the non-spectator contributions may be important JJ for 
evaluating the lifetimes of B and D mesons. However, following literature, in the exclu- 
sive processes, the spectator seems to be fully dominant and it is reasonable to believe 
that the evaluation is sufficiently accurate if only the spectator contribution is taken into 
account. 

In the quark-diquark picture, the two-body baryonic wavefunction can be written as 


i)i(x x , kj_) = A^xiX2exp[-6 2 (k^_ + mj(xi - x 0i ) 2 )}. (5) 
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where i = A c , n,p or A, x±, x 2 {x 2 = 1 — 27) are the longitudinal momentum components 
of the quark and diqaurk, kj_ is the transverse momentum of the constituents, and TV* is 
a normalization factor. 



The two form factors F\ and F 2 in eq.(|j) are related to an overlapping integral of the 
wavefunctions of A c and the produced baryon ||13|| : 



where /(k>) is the overlapping integral of the wavefunctions of A c and / 
m = f^-) 7/ %- 9/2 [A / ir 6 (v / 26 e )]- 1/2 expf-26 2 e ^" 1 



u + lj /J V u+ 1 

^-Jtt . , „/ 2be \ ( 2be 



2bi 



1 dz exp( -^ 2) I y - t^tt " "J r + ^TlJ ' (7) 



where y = bmf\Juj + 1, is the velocity transfer and u — v ■ pf/rrif. Af and are 
defined as 

rl 
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Af = I dx x 6 (l — x) 2 exp[— 2b 2 m 2 (x — e/rrif) 2 ], 

j 

K 6 (V2be) = / da; exp(— x 2 )(x + 



-yfibt 

It is noted that here the limit m c — > 00 is taken. In general, according to the Lorentz 
structure the transition matrix element of -Bj(pi) — > Bf(jpf) is decomposed into the form 

-(^i(? 2 )7m + m(q 2 )cr^q u + g-z{q 2 )q^}Ui(pi), (9) 
where q = Pi — Pf- Comparing eqs. (01) and (|9]) we have relations 

/i(9 2 ) = 9i(q 2 ) = F 1 (q 2 ) + ^F 2 (q 2 ), (10) 

rrii 

/ 2 (g 2 ) = 92{q 2 ) = h{q 2 ) = 9,{q 2 ) = —F 2 (q 2 ). (11) 
Finally we reach the following formula: 

{B f {p f )\Vp - A^Biipi)) = U f [ax^ 5 + a 2 { Pi )^ + &i7a« + hipd^Ui, (12) 



where 

ai = -[Fi(g 2 ) + F 2 {q% a 2 = -2F 2 {q 2 )/ mi ; 

h = [F^q 2 ) - F 2 (q% b 2 = 2F 2 {q 2 )/m i . (13) 



The process A c — > pn° is Cabibbo-suppressed and the reaction corresponds to the 



so-called inner emission which is a color-suppressed process |TlJ. In terms of the weak 
effective hamiltonian we obtain 



M(A C 



pir 



-^v ub v; b 



1 + — )(c 3 + c 4 )(7r°|^(l - 7b)«|0)<p|«7m(1 - 7 5 )c|A c > 



+ (c 5 + ^)(vr°|^(l + 75 )«|0)(p|w7 M (l - 7s)c|A c ) 
-2(^ + c 6 )(tt |«(1 + l5 )u\0)(p\u(l - 75 )c|A c ) 

+ ^K^( Cl + ^)(vr°|d7 M (l - 75)rf|0}(p|«7M(l - 7 5 )c|A c ) 



Gf - 



-^v ub v: b 



(c 3 + ^)(tt |J 7 m (1 - 75 )rf|0)(p|^ 7M (l - 75 )c|A c ) 



+ (c 5 + ^)(vr°|J 7 M (l + 7fe)d|0>(p|«7 M (l - 7s)c|A c ) 



(14) 



where N c is the effective number of color which includes non-factorizable effects. After 
some simple manipulations, it is recast as 



M(A C 



PIT 



Gf 
V2 



V ud V: d ( Cl + + V ub V* b (^j- + c 4 )] ^U P (A> + B' l5 )U Ac 



C3_ 



+ V2G F V ub V; b (^- + ce)(-^|) U P (A" + B" l5 )U Ac 
tU p (A + B l5 )U AcJ 



if* 



(15) 



where 
A = 
B = 
A" 



Gi 



^-U{ [v ud v; d ( Cl + ^) + v ub v^ + c 4 )]a' + v ub v; b (^- + c 6 ) a"} 



c 3 



G, 



c 2 



N c 
c 3 



c 5 



'A", 



c 5 



[K^(ci + ^) + K 6 v c * 6 (^- + C4 )]b> + v ub v; b (^ + C6 ) b"}, 



mz 



m u (m c - m u ) 



A' , B" 



mt 



m u (m c + m u ) ' 
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m, 



+ mz — m 



A' = bi(m Ac - m p ) + b 2 r> 
ai = -[Ftiml) + F 2 {ml)}, a 2 = -2F 2 {ml)/m A 
h = [F^ml) - F 2 {ml)}, b 2 = 2F 2 (m 2 n )/m Ac . 



P -,B> = - ai (m Ac +m p ) 



m. 



+ mi — mi 



a 2 - 



(16) 



The expressions for A c — > mr + (p + ) are similar to that for A c — > pir , so we omit the 
details and show the result directly, 



M(A c ^nn + ) = iU n {A + B lb )U A 



(17) 



here A, B take values as follows 

A : 



+ 02) - + C4)]AA' 



A', 



5 



N c (m u + m d )(m c - m d ) 



V2G F f n m 2 n 



B\ 



+v ub v; b (^ + c 6 y 

N c (m u + m d ) (m c + m d ) 

where A' and B' have the same forms as those in eq.([l(J) with m p being replaced by m n . 
In the same way, we have 



M(A C -> np + ) = 



- VUVJt^- + C4)](/|K7"rf|0)<n|<i7„(l - lMA c ). (19) 



Finally, 



A4(A C -> np H 



Z7 e* M 



^i7 M 75 + ^2(paJ m 75 + ^i7 M + B 2 (p Ac )i 



where Ai = aii], Bi = bm (i=l, 2), and 

Gp 
V2 



V = + c 2 ) - V ub V; b (§- + c A )]f p m p . 



'Nr. 



Ua c , (20) 



(21) 
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For A c — > AK + (K + *), the formulas are in analog to those given above. We have 



M(A C — > AK + ) 
V2 



c 3 



'iV, 



'-/V, 



c 5 



iV2G F f K m 2 K 



(m u + m s )(m 2 c - m 2 ) 



U A (m c + m s )A' - (m c - m s )B'^y 5 U A 



= iU A (A + B l5 )U Ac , 
where A, B are given by 



(22) 



a = ^[v us v; s (^ + c 2 )-v ub v; b (^ + CA )}f K A', 

\r t/v Cs j_ a V2G F f K m 2 K 

-VubV cb (— + c 6 - ■ - -A , 

N c (m u + m s )(m c - m s ) 



B 



(m u + m s )(m c - m s ) 
= ^[Vu S V;^ + c 2 )-V ub V: b (^- + c 4 )]f K B' 

4-V v*(° 5 | „ \ V2G F f K m 2 K 

+ v ubV cb {— + c 6 J- ■ ■ -B , 

N c [m u + m s ){m c + m s ) 



(23) 



with 



a ,,,2 n/ -2 



m A c + m K 



A' = b 1 (m Ac - m A ) + b 2 (p K • PaJ = 6i( m A c - m A ) + b 2 

B = a 1 (-m Ac - m A ) + a 2 (p K ■ p A J = -ai(m Ac + m A ) + a 2 c - .(24) 

For A c -> Aif" + , we have 

M(A C ^A1T + ) = U A e*^\A a ^ 5 + A 2 (p A )^ 5 + B 1 ^ + B 2 (p A ) lx ]u Ac , (25) 



where A* = and B$ = bifj (i=l, 2), and 



77 = -£[K,V^(§- + c 2 ) - K&KUtt + c 4 )]/^m^. 



(26) 



With all the amplitudes, we can immediately obtain the decay widths and correspond- 
ing branching ratios if there were no final state interaction. But in fact, the FSI effects 
can drastically change the results. In the following, we will present the formulas for the 
FSI effects. 
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Fig.l The Final State Interaction for A c — > pn 1 



Since the final state interactions may change the whole picture for some processes, 
one needs to take into account these effects seriously. Unfortunately, the FSI processes 
are governed by non-perturbative QCD, so we cannot simply start from any underlying 
theory to carry out the computations yet. One could extract valuable information from 
data. Theoretically, we may also use some phenomenological approaches to study the 
order of magnitude of the FSI effects. 

There are several ways to evaluate the FSI. A usually adopted method is the one- 



particle exchange model |T4| and another one is that in which the produced primary 
particles of the decay scatter into the final products via exchanging the Regge poles |T!| . 
The first method is straightforward and simple, and we employ it in this work. 

In this method, the effective vertices can be obtained from the chiral Lagrangian |16| . 
The coefficients in the chiral Lagrangian are obtained directly from data where all external 
particles are on their mass shells. Instead, for our case the exchanged mesons or baryons 
are obviously off mass shell. To compensate the off-shell effects, one needs to introduce a 



form factor at the effective vertices 17 



The Feynman diagrams for the FSI are shown in Fig.l. 
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In fact, there are some u-channel diagrams, namely in Fig.l, p and tc° may exchange 
positions and the intermediate bosons p(K*) or ir(K) are replaced by baryons p or A. But a 
direct calculation indicates that the u-channel contributions are only at most 1CT 3 ~ 10~ 2 
of the t-channel, so we omit the u-channel part in later calculations for simplicity. 

Moreover, the direct decay of A c — > p7i° is color suppressed compared to A c — > nn + etc., 
and hence has smaller amplitudes. Thus we only consider A c — ► tot + , np + , AK + , AK + * — > 
pn , but not the opposite directions. 

For A c — > pn , the amplitude with final state interactions can be written as 

AFSI 1 f d3 Pl d3 P? fn \^( , \ 

■A{A c ^nn + )(nn + \S\ P 7r°), (27) 

where for clarity, we set p\ = p n ,P2 = Pn+,P3 = Vy an d P4 = Pit - In terms of the Cutkosky 
cutting rule 0|, we can easily calculate the inelastic scattering amplitude. 



The effective pNN has the following form [E2 



£ p nn = QpNN^N^ + Kp ^^)N. (28) 



With this vertex, we have 



< a (p2+PA) a (p'i + rm)i(A + B l5 )U Ac ( PAc ). (29) 



Here F(k 2 ) is a form factor |20|, |2l |. In fact, the effective couplings at the vertices are 
obtained from real processes where all particles are on mass shell, but in our case, at 
least one of the particles (the exchanged p meson) is off shell. Therefore, we introduce 
this form factor as a compensation to the off shell effects. Following the literature, in our 
calculation we take 

F(k 2 ) = [(a 2 -m 2 p )/(a 2 -p 2 p )}\ 

where a is a parameter. 



10 



The explicit expressions for A^+_^ as well as those for other inelastic final state 
interactions, A c — > np + — > p7r°, A c — > AK + — > p7r° and A c — > AK* + — > £>7r°, are rather 
complicated, so we collect them in the appendix. 

III. The numerical results 

As inputs, we take the q's in the weak effective Lagrangian from literature |], [Tj| where 
the renormalization is carried out to one-loop level: 

Ci(m c ) < Oi(m c ) >=c[< Oi > tree } 

where d i are effective Wilson coefficients. 

In the wavefunction ipi(xi, kj_), we take b = 1.77 GeV and 1.18 GeV corresponding to 
< k 2 ± > l / 2 = 400 MeV and 600 MeV, respectively. x i{x 0i = 1 — m D /mi, i = A c ,n or A) 
is calculated with the diquark mass mu being taken the value of 600 MeV. We let the 
parameter a in the form factor F(q 2 , a) vary in a reasonable region. The results are not 
very sensitive to these parameters as shown below. 

We also set m u ~ ~ 10 MeV, m s = 150 MeV, m c = 1.35 GeV. For the hadron 
decay constants, f n = 132 MeV, f K = 160 MeV, f p = 210 MeV, and f K * = 221 MeV. 

The CKM entries are in the Wolfenstein parametrization and all values are taken from 
the Particle Data book. 

a. Without the final state interaction. 

We first present the numerical results when the FSI effects are turned off. The numer- 
ical results are shown in Table. 1 and Table. 2. 



Products 


P7T° 




np + 


AK + 


AK* + 


b= 1.77 

GeV 




0.207 


0.209 


0.264 


0.301 


0.384 


F 2 


-0.0533 


-0.0538 


-0.0730 


-0.0728 


-0.0989 


b= 1.18 
GeV 


Ft 


0.145 


0.146 


0.165 


0.208 


0.236 


F 2 


-0.0372 


-0.0375 


-0.0456 


-0.0503 


-0.0609 



Table 1. The form factors F\ and F 2 as b takes two values. 
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Branching ratos 


p7T° 


U7l + 


np + 


AK + 


AK* + 


b = 1.77 GeV 


4.35 x 10" 6 


2.10 x 10" 4 


3.57 x 10~ 4 


3.86 x 10" 4 


4.91 x 10" 4 


b= 1.18 GeV 


2.14 x 1(T 6 


1.03 x 10" 4 


1.39 x 10~ 4 


1.84 x 10" 4 


1.85 x 10" 4 



Table 2. The branching ratios. Here we have ta c = 0.206 ps [27|. 



b. Taking the FSI effects into account. 

The coupling constants of the strong interaction at the vertices take the values gK*Kn = 
5.8, g pnw = 6.1 fH g6f, g pNN = 2.5, n p = 8.0 @, g vNN = 13.6 @, g KAN = -10 0. 
To determine the coupling gK*AN, we assume the equality gx*AN I 'g P NN = gKAN / g-wN n ■ 
In the compensation form factors F(a), the parameter a takes values in the reasonable 



region]!?], 20), Q as 1.2 — 2.0 GeV. In our concrete calculations, we employ a = 0.9 GeV, 



1.0 GeV, and 1.2 GeV, respectively. From we have \gKAN / '9khn\ = 12, thus we do 
not need to consider the contribution of A c — > UK to A c — > pir via final state interaction. 
The calculations are done for the case of maximum interference. 
The numerical results with FSI effects are listed in table 3. 





without 
FSI effects 


with FSI effects 


a = 0.9 GeV 


a = 1.0 GeV 


a = 1.2 GeV 


b = 1.77 GeV 


4.35 x 10~ 6 


3.14 x 10~ 4 


2.92 x 10" 4 


3.55 x 10~ 4 


b = 1.18 GeV 


2.14 x 10~ 6 


1.20 x 10" 4 


1.11 x 10" 4 


1.50 x 10" 4 



Table 3. The branching ratio of A c — > p?r with and without the FSI effects. 



IV. Conclusion and discussion 

It is widely recognized that non-perturbative QCD is, so far, an unsolved problem, but 
it is entangled with the effects of the fundamental mechanisms. We need to explore its 
properties in the phenomenological studies and gradually get better understanding of it. 

The processes where only mesons are involved have been carefully investigated, includ- 
ing the final state interactions, even though there is still no way to evaluate them exactly 
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from underlying theory. On the other hand, for the baryon case, the situation is more 
complicated because there are three valence quarks in a baryon. 

It is reasonable to consider the diquark structure in baryons, which may be one of the 
possible configurations of three quarks. For a heavy baryon with one heavy quark and 
two light quarks, the two light flavors constitute a color 3 diquark where the interaction 
is attractive, and the heavy one acts as a color field source. This configuration may be 
expected as the dominant one. Thus in this work we employ this physical picture to carry 
out the calculations. Definitely, this configuration does not possess 100% probability, so 
we might introduce a phenomenological parameter which is less than unity, but close to, 
for describing the deviation. Therefore, the branching ratios we obtained in the text must 
be multiplied by this parameter. However, the relative ratios of A c — > pir° to A c — > niT + 
does not depend on it and this fact can be tested in experiments. As a matter of fact, 
some authors also used the diqaurk picture to evaluate the decays of heavy flavors to 
lighter ones recently f25|j . 

Moreover, it is easy to accept that in the spectator picture, the heavy quark turns into 
a lighter one, there could be a possibility that a light quark in the light diquark breaks 
out and combines with the newly produced quark to constitute another diquark. But 
from the physics intuition, this probability should be smaller than retaining the original 
diquark structure as spectating the heavy quark turn into a light one, so we assume that 
the light diquark serves as a spectator and does not change after the weak transition. 

When we calculate weak transitions, we use the factorization ansatz, which is not 
accurate. Therefore the absolute values of the obtained branching ratios may decline 
from the real values by a certain amount. But according to the information we have 
learned in the meson case, for such exclusive processes, factorization works. We will 
pursue this problem in our later works. 

In this work, we show that the FSI effects are important, not only for the meson cases, 
but also for baryon-involved processes. It is indicated that without considering the FSI, 
the estimated branching ratio of A c — > pir° is almost two orders smaller than that of 
A c — > mr + , but with FSI, they could have the same order. The numerical results show 
that the theoretically calculated values depend on the parameters a and b in the model 
employed here. This originates from the lack of knowledge of non-perturbative QCD, and 
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we cannot fix them at present. However, as long as they fall in a reasonable region, the 
order of magnitude does not change. 

For evaluating the FSI effects, we only consider the absorptive part of the triangle 
diagram. It is the real final state process because both hadron which are produced via 
the weak transition are on their mass shell. But definitely, the dispersive part of the 
triangle which needs to be renormalized would also bring up some effects and we ignore 
them in this calculation. Higher order FSI may also play roles and it can induce a strong 
phase and cause CP violations. We have investigated the possible CP violation in the 
meson case p0| , |21|| because there are data about the scattering amplitude and phases. In 
this concerned case, it is impossible to carry out accurate calculation yet, so we can only 
trust it to the order of magnitude. 

By including the final state interaction, we show that the branching ratios of A c — > pn 
could be greatly enhanced by almost two orders and it is similar to the meson case, as 
D° — > K°K° is much suppressed compared to D° — > K + K~ , but measurements indicate 
that they are at the same order. This is obviously due to the FSI effects. Therefore when 
we evaluate the processes for baryons, we cannot ignore the FSI, as we understand in the 
meson case. 

As pointed out above, many uncertainties exist in the theoretical calculations, so for 
a better understanding, we need more data as well as a more reliable theory for the non- 
perturbative QCD. Fortunately, the BEPC, CLEO and other high energy accelerators 
such as the B-factory and even TEVATRON and LHC will provide much more data in 
the future and we will be able to use them to determine the parameters involved in the 
phenomenological models and surely test the model itself. 
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Appendix 

In this appendix we present some details of the formulation about the final state in- 
elastic scattering amplitudes. 
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(1) For A c 
A? 31 



117T^ — > p7r , 

d4pi -5(pl - m 2 )5{p\ - m ^ k2 )9 P NN9^ 



k 2 



- 4pa c • Pi + m\ -{m\fa - m\i>K c ~ %Pa c • Pi fa) 

m N 

+mi 



m i( m 2 — m V' 



(2pA c - ps) + — ^-(P2 • Pa — M 2 + 2p Ac ■ Pl + fa fa 
m N 

{m22 ~ ml) -fa]^(A + B l5 )U Ac (PA c ) 



1 f d 4 Pl 
2 



8{p\ - m 2 )5(p\ - ^^^SpNNg^ 



k 2 — m 2 



(2tt) 2 

U p {p 3 )i(AH? + BH^ 5 )U Ac ( PAc ); 

1 dMWTO,^ + BH a j5)UAc{pAc 



i 16ttM 



A; 2 



(30) 



where 9 is the angle spanned between pi and P3, k 2 = (pi — p 3 ) 2 = m\ + m 2 — 2E\E 3 + 
2| p\|| p" 3 | cos(6>), and \~f x \ = A 1/2 (M 2 , m 2 , m|)/(2M). A(a;,y,z) = x 2 + + z 2 - 2xy - 
2?/2 — 2a;^ is the well-known function. In the above expression, A, B take the values in 
Eq.flTSD and iff, #f are 



- 4p Ac • pi + m 1 



m N 



■{m\m 3 - m\M - 2p Ac ■ pim 3 ) - 



m l( m 2 — m 1) 



+mi 



2M + m 3 + — ^(p 2 -P4-M 2 + 2p Ac ■ p x + m 3 M) + 



(m 2 , — m 2 )m 3 



m. 



H a 2 



4pA c • Pi + m? — — (m?m 3 + m?M - 2p Ac ■ Pim 3 ) 



m 2 (m 2 — m 2 ) 



+mi 
(2) For A c - 



rrip 

,2 



2M + m 3 + -^-(p 2 ■ P4 — M 2 + 2p Ac • px - m 3 M) + 



(m 2 , — 1711)1713 



171, 



pn we have similar expressions 



np+— >p-7r° 



oPpx d 3 p 2 



2 J (2tt) 3 2E 1 (2tt) 3 2E 2 

5{p\ — m 2 )5{p\ — m 2 



d^pi 



(2vr)V(p 1 +p 2 -p Ac ) ■ A(A C - nn + )(np + \S\pn°) 
iF{k 2 )g 

71 



(2tt) s 



fc 2 



-iU p (p 3 )j 5 (fa + mi) 



• e p( fc + P4)/" [^17^75 + A 2 (pi) v 7 5 + S l7 , + ^(pO.Jc/a. 
where p\ = p n ,p 2 = P P +,P3 = P P and p 4 = p n o. Eq.(|31[) is recast as 



(31) 



4 

np+— ►p7r° 



'^(p 2 - m 2 )5(p? - m 2 ) ^ (fc2 2 } ^f^ ^ p (p 3 )[A + B bl5 }U Ac 

k z — mi 
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1 d{cosd)\~f 1 \iF{k 2 )g^ NN g fm 



16nM 



k 2 — m 2 



~-iU p (p 3 )[A b + B bl5 ]U Ac , 



(32) 



and A b , B b are defined as 

A b = 2m 1 {A 1 (M - m 3 ) + A 2Pl ■ (p 3 - p Ac ) + 

P2-P4 2 



P2 'Pi 



rrin 



-A 1 M + A 2 pi ■ (pAc-Pi)]} 



+2{2 Pl -p 3 A t 



ml 



-mlA,}, 



B h = 2m l {B l {-M-m 3 )+B 2 p 1 -(p 3 -p Ac ) + ^^\B 1 M + B 2 p 1 -{ PAo -p 1 )]} 



ml 



+2{2p l -p 3 B 1 - 



P2-P4 2 



ml 



mlB,}. 



and At = a t m B { = (i = 1,2). 

(3) For A c — > A7a" + — > p7r°, we have similar expressions 



a fsi 

1 r d 3 pi d 3 p 2 

2 J (2ir) 3 2E 1 (2vr) 3 2£ 2 
1 d(cos6>)| p\| iF(k 2 )g K * NA g K *+ KlT - 



-i 16vrM 

where A, I? take the values in Eq.(p3|) and if^, if| are 



(2vr) 4 ( 5 4 (p 1 + p 2 -pa„) • ^(A c - Air* + )(AK + |5|j97r°) 

^(pa)^^ + BHfr 5 )U Ac (p Ac ), (33) 



A; 2 — m 2 ^- 



A? 



HI 



- 4 Pa c •pi+m 1 - 
+m l \-2M + m 3 + 

- 4p Ac • Pi + ml - 
+m l \2M + m 3 + 



2k' 



m N + m A 

K 

m N + m A 

K 

m N + m A 
2k' 



{m\m 3 - m\M - 2p Ac • Piw 3 ) - 



m l( m 2 — m 4) 



(p 2 ■ p 4 - M 2 + 2p Ac • p! + m 3 M) + 



(m 2 — m|)m 3 - 



(m\m 3 + m 2 M - 2p Ac • pim 3 ) 



mf (mj — m|) 



A'* 



m^v + m A 



(p 2 • p4 — M + 2p Ac • p! - m 3 M) + 



(wij — m 2 )m 3 



m 



A'* 



where «' is from the interaction vertex of fl28|) and we have pi = p A ,p 2 = Pk*+,P3 = P P , 

P4 = Ptt°- 



(4) For A c -> A7T+ -> pvr°, 

/1 AA'*+^p7rO 

1 /■ c? 3 pi d 3 p 2 



2 y (2vr) 3 2£i (2tt) 3 2E 2 



(27r) 4 5 4 (p! +p 2 -p Ac ) • A(A C - AiT + )(Air+|,S|p7r > 



1 d(cogg)[gjJ iF(k 2 )g KAN g K , 
-i 16vrM A; 2 - 



^^3)^ + ^75]^, (34) 
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where p\ = p A ,p 2 = Pk*+,P3 = P P and p 4 = p n o. A d , B d are defined as 



A d = 2m 1 {A 1 (M -m 3 )+A 2 p 1 -( P3 - PA J + ^p\-A 1 M + A 2 p 1 -(p Ac -p 1 )]} 
1 m 2 l 



+2{2 Pl ■ p 3 A } 



P2-P4 2 



ml 



B d = 2m 1 {B 1 (-M-m 3 ) + B 2 p 1 -(p 3 -p Ac ) + ^p\B 1 M + B 2 p 1 -(p Ac -p 1 )]} 
+2{2 Pl ■p^-^pmjBA, 



mij 



and A, = a^, ^ = brf (i = 1,2). 
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